Level statistics in a metallic sample: Corrections to the Wigner— Dyson distribution. 



Vladimir E.Kravtsov^f and Alexander D. Mirlin^J 
^ International Centre for Theoretical Physics, P.O. Box 586, 34100 Trieste, Italy 
^ Institut fiir Theorie der Kondensierten Materie, Universitdt Karlsruhe, 76128, Karlsruhe, Germany 

(June 23, 1994) 



ON 

^ ■ 
m 



> 

O 

00 

O 

ON 



o 



Deviation of the level correlation function in a mesoscopic 
metallic sample from the Wigner-Dyson distribution is calcu- 
lated by using a combination of the renormalization group and 
non-perturbative treatment. For given spatial dimension the 
found correction is determined by the sample conductance. 
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The problem of level correlation in quantum systems 
attracts interest of physicists since the work of Wigner 
The random matrix theory developed by Wigner 
and Dyson [|| describes well the level statistics of various 
complex systems, like nuclei or atoms. Later on, Gor'kov 
and Eliashberg put forward an assumption that the 
random matrix theory is also applicable to the problem 
of energy level correlations of a quantum particle in a 
random potential. This hypothesis was proven by Efetov, 
who showed Q| that the level-level correlation function 
R{u}) (its formal definition is given below) is described 
by the Wigner-Dyson distribution for uj <^ Ec, Ec being 
the Thouless energy. For uj ^ Ec, the behavior of the 
correlation function changes because the corresponding 
time scale lj"^ is smaller than the time E~^ the particle 
needs to diffuse through the sample. The form of the 
correlation function in this region is dependent on spatial 
dimensionality and was calculated in Ref. by means 
of the diffusion perturbation theory. 

In the present Letter we find a correction to the 
Wigner-Dyson distribution in the region w ~ A <C -Be, 
A being the mean level spacing. This is not a trivial task, 
because we calculate a correction to the result which is 
essentially non-perturbative. 

We study the two-level correlation function R{s) de- 
fined as 



i?(.) = 



1 



{i^{E)i^{E + u)) 



(1) 



where s = ui/A, v{E) is the density of states at the en- 
ergy E and (. . .) denote averaging over realizations of the 
random potential. As was shown by Efetov the corre- 
lator (1) can be expressed in terms of a Green function of 
certain supermatrix cr-model. Depending on whether the 
time reversal and spin rotation symmetries are broken or 
not, one of three different cr-models is relevant, with uni- 
tary, orthogonal or symplectic symmetry group. We will 
consider the case of the unitary symmetry (correspond- 
ing to the broken time reversal invariance) throughout 
the paper; the results for two other cases will be pre- 



sented at the end. The expression for R{s) in terms of 
the (T-model then reads: 
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Re / DQ{r) 



expj-^ j d'^rSiT:[-D{VQf -2iLokQ]^ 



(2) 



Here Q = T^^AT is 4 x 4 supermatrix, with T belong- 
ing to the coset space U{1, 1|2), A = diag{\, 1, —1, —1}, 
k — diagjl, — 1, 1, — 1}, Str denotes the supertrace, V 
is the system volume and D is the classical diffusion 
constant. To find the detailed description of the model 
and mathematical entities involved, a reader is refered to 
Refs. |§. 

For w -C i?c ~ D/ (L being the system size, so that 
V = L"^) the leading contribution to the integral (||) is 
given by the spatially uniform fields Q{r) — Q. Then the 
functional integral in eq. (|^) reduces to an integral over a 
single supermatrix Q and can be calculated yielding the 
Wigner-Dyson distribution M: 



Rwd{s) = 1 - 



sin^(7rs) 



(3) 



The aim of the present paper is to calculate a correc- 
tion to eq.(||) due to spatial fluctuations of Q{r) in 
eq.(||). The procedure we are using for this purpose 
is as follows. We first decompose Q into the constant 
part Qo and the contribution Q of higher modes with 
non-zero momenta. Then we use the renormalization 
group ideas and integrate out all fast modes. This can 
be done perturbatively provided the dimensionless con- 
ductance g — Ec/A :s> I. As a result, we get an in- 
tegral over the matrix Qq only, which has to be calcu- 
lated non-perturbatively. We believe that this combina- 
tion of the perturbative renormalization-group-type and 
non-perturbative treatment is of a methodological inter- 
est and might be used for other applications. 

To begin with, we present the correlator R{s) in the 
form 



R{s 



1 



J DQexp{-J^{Q)}\u^o 



(27ri)2 du^ 

T{Q) = -\ j Str(VQ)2 + 5 j StrAQ + u j StrQAfc (4) 



where 1/t = ttvD/A, s — ixsjliV , u - 
decompose Q in the following way 

g(r) = To-ig(r-)To 



■Ku/2iV . Now we 



(5) 
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where Tq is a spatially uniform matrix and Q describes 
all modes with non-zero momenta. When uj <^ Ec, the 
matrix Q fluctuates only weakly near the origin A of the 
coset space. In the leading order, Q = A, thus reducing 
to a zero-dimensional cr-model, which leads to the 
Wigner-Dyson distribution (^). To find the corrections, 
we should expand the matrix Q around the origin A. This 
expansion starts as Q 

= A ( 1 + T4^ + — + . . . I , (6) 



where is a supermatrix with the following block struc- 
ture: 



W 



ti2 
t21 



(7) 



Substituting this expansion into eq, 
T = To + Ti + 0{W^) ; 

= / Str 



we get 



-{VWf + SQoA + uQQAk 



= \ j Str[SC/oATy2 + nUokAW"^ 



(8) 



where Qo - T^^AT^, f/o = ToKT^\ U^k = ToAkT^\ 
Let us define ^e//(Qo) as a result of elimination of the 
fast modes: 



g-^e//(Qo) ^ g-.^o(Qo)^g-.^i(Qo,H^)^^ 



(9) 

where (. . .)w denote the inetgration over W. Expanding 
up to the order W^, we get 



(10) 



The integral over the fast modes can be calculated now 
using the Wick theorem and the contraction rules |^,^ : 

{StTW{r)PW{r')R) = n(r - r') 

x(StrPStri? - StrPAStri?A) ; 
{Str[W{r)P]Sti[W{r')R]) = n(r - r')Sti{PR - PARA) 

nir) ^ I (11) 



Using now eq. (|4|) , we get the following expression for the 
correlator to the order: 



X[l- 

8g2 



16,g2 



s"(StrQoA)^ 



(StrgoAA:))^ 



Ts(StrQoA)(StrQoAfc)^ 



X exp{— — sStrQoA} 
2i 



(14) 



This integral over the supermatrix Qo can be calculated 
by using the known technique 0] , yielding 



R{s) = 1 - 



sin^(7rs) ad 
(7rs)2 ^ 



sin^(7rs) 



(15) 



The last term in eq. ([l5|) just represents the correction 
of order 1/g^ to the Wigner-Dyson distribution. The 
formula (^_5|) is valid for s <^ g. Let us note that the 
smooth (non-oscillating) part of this correction in the 
region 1 ^ s ^ g can be found by using purely pertur- 
bative approach |^,^ . For s ^ 1 the leading perturbative 
contribution to R{s) is given by a two-diffuson diagram: 



R 



pert (s) - 1 = 2^ RC 

qi = Tirii/L 
n,^ 0,1,2,.. 

1 

Hp > 

27r2 



1 



ni>0 



(—is + TT^gn^) 



(L>g2 - iujf 

(16) 



At s ^ 5 this expression is dominated by the q = term, 
with other terms giving a correction of order 



R. 



■pert 



(s) = i 



1 



ad 



27r2s2 



27r252 



(17) 



{2t:Y TTvDq^ 

where P and R are arbitrary supermatrices. The result 
is: 

(^i) = ; 

i^i)^^ J drdr'n^ir -r'){:sStYQoA + uStrQoAk)^ (12) 



Substituting eq.(12) into eq.(|lCl|), we find 



J^effiQo) = — sStrQoA+ — uStrQoAfc 



where was defined in eq. (y^ . This formula is obtained 
in the region I <^ s ^ g and is perturbative in both 1/s 
and 1/g. It does not contain oscillations (which cannot 
be found perturbatively) and gives no information about 
actual small-s behavior of R{s). The result (^5|) of the 
present Letter is much stronger: it represents the exact 
(non-perturbative in 1/s) form of the correction in the 
whole region s <^ g. 

The calculation presented above can be straightfor- 
wardly generalized to the other symmetry cases. The 
result can be presented in a form valid for all three cases: 



165= 



(sStrQoA -1- -uStrQoAA:)^ ; 



Rp{s) 
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(13) 



where /3 = 1(2,4) for the case of orthogonal (unitary, 
symplectic) symmetry; R^^ denotes the corresponding 
Wigner-Dyson distribution, explicit form of which can 
be found e.g. in HJ^. ( We denote by g the conductance 



2 



per one spin projection: g = Ef-j IS. = vDL'^~^, without 
multiphcation by factor 2 due to the spin.) 

For s ^ the Wigner-Dyson distribution has the fol- 
lowing behavior: 
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(19) 



As is clear from eg. ([iq) , the found correction does not 
change the power /3, but renormalizes the prefactor C/j: 



Rp{s) = 1 



(/3 + 2)(/3+l) ag 
2/3 TT^g"^ 







(20) 



As was already noted, the above derivation of eq.dl^) 
is valid provided g 3> 1, i.e. in the case of a good metal. 
The level statistics is also of big interest in the criti- 
cal point g = 5* corresponding to the Anderson metal- 
insulator transition. The "tail" of this distribution at 
s ^ g was studied in the recent paper ||lo[] , but nothing 
has been known analytically about its behavior at s ^ 5. 
In c? = 3 the critical point corresponds to the strong 
disorder: 17* ~ 1, where the Wigner-Dyson distribution 
loses its validity. However, we can consider formally the 
spatial dimension to be d = 2 -I- e with e <C 1, where 
(7* ~ 1/e 3> 1, i.e. the transition point is still in the weak 
disorder regime, (one may hope that the obtained results 
will be still qualitatively valid when continued to e=l.) 
Then the Wigner-Dyson distribution has again a para- 
metrically broad range of validity s < 5* , and a question 
about the corrections to it due to spatial fluctuations is 
meaningful. We can repeat basically the same calcula- 
tion, as was done above for the d — 3 metallic system. 
However, the diffusion propagator n(r) should now be 
replaced by its form in the critical region ||Tll|: 



H(r) oc 



d'^q exp(zqr) 



(2n) 



9*1 



As a result, we find in close analogy with eq.(15), 



1 



SR{s) oc —J- sin (tts 



(21) 



(22) 



It should be noted that in contrast to eq.(p^, the ex- 
act numerical coefficient in eqs.(^l|), (|2^ ) is not known. 
Eq.(p2[) is written for the unitary symmetry case P — 2; 
the general expression for the correction is given by 
eq.dlq), with substituted for g. 

As is seen from eqs.(|20|), (|22^, the found corrections 
lead to an enhancement of the prefactor cp in the low- 
s behavior R{s) cps^ of the correlator, that means 
physically a weakening of the level repulsion. The cor- 
rection to the prefactor cp is of order 1/g^. When a sys- 
tem approaches the metal-insulator transition in d = 3, 
the conductance g tends to a value 5,, which is of order 



unity. The calculation carried out in the present paper 
is no more justified in this region. We expect, however, 
that being formally applied it gives qualitatively correct 
behavior of R{s) at small values of s, namely an enhance- 
ment of the coefficient cp by a factor of order unity. 

In conclusion, we have calculated the deviation of the 
level-level correlation function Rp{s) in a mesoscopic 
metallic sample from its universal Wigner-Dyson form, 
using the supersymmetric sigma-model approach. The 
method of calculation combines a perturbative elimina- 
tion of fast spatial modes (in spirit of renormalization 
group ideas) and subsequent non-perturbative evaluation 
of an integral over the zero mode. The found correction 
is of order l/g^, where g is the dimensionless conduc- 
tance. It does not change the power /3 of the behavior 
of the correlator R{s) as s ^ 0, but renormalizes the 
corresponding prefactor. 
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